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A  newly  developed  theory  for  predicting  the  response  of  a  linear  system  to  an 
electromagnetic  pulse,  based  only  on  the  measured  continuous-wave  magnitude, 
is  applied  to  the  problem  of  possible  electromagnetic  interferences  at  a  sensitive 
part  of  a  torpedo.  The  measured  magnitude  representing  the  system's  transfer 
function  is  deduced  first  from  the  measured  response  at  this  sensitive  point  to  a 
known  cw  source,  supplied  by  the  Naval  Surface  Warfare  Center.  We  derive  an 
analytic  expression  for  the  magnitude  square  of  the  transfer  function  to 
approximate  the  measured  data  and  obtain  a  system  transfer  function  in  terms  of 
the  complex  fi"equency,  from  which  we  predict  the  system's  cw  phase 
characteristics  and  its  multiple  solutions  due  to  a  given  impulse  source. 

Key  words:  convolution  integral;  impulse  response;  linear  system;  magnitude; 
phase;  transfer  function. 


1.  Introduction 

A  theoretical  method  was  developed  to  predict  the  frequency  and  time  characteristics  of  an 
unknown  linear  system,  based  only  on  the  measured  cw  magnitude  data,  without  requiring  phase 
information  [1].  The  approach  is  to  deduce  an  approximate  analytical  expression  for  the  overall 
system  magnitude  square  from  the  given  measured  magnitude  data.  The  essential  steps  involved 
in  achieving  this  purpose  [1]  may  be  summarized  as  follows:  (i)  identify  the  resonant  frequencies 
and  the  maximum  responses  at  these  frequencies  from  the  data,  (ii)  approximate  each  of  the 
resonant  regions  by  a  component  magnitude  square  with  a  simple  second-order  transfer  function, 
(iii)  sum  up  all  of  the  component  magnitude  squares  to  yield  an  overall  magnitude  square  to 
represent  the  unknown  system,  (iv)  derive  from  it  a  set  of  multiple  transfer  functions  in  terms  of 
the  complex  frequency,  (v)  determine  the  corresponding  phases  and  impulse  functions,  and  (vi) 
compute  the  time  responses  of  the  system  to  a  known  electromagnetic  pulse  (emp)  source  by 
performing  analytically  the  convolution  integral. 

Steps  (i),  (ii),  and  (iii)  are  outlined  in  section  4.1.  The  approximation  for  each  resonant 
frequency  by  a  second-order  transfer  function  as  required  in  step  (ii)  is  presented  in  section  2.  A 
short  discussion  on  frequency  transformation  for  avoiding  manipulations  with  large  numbers  is 
given  in  section  3.  The  derivation  of  transfer  functions  from  the  approximate  magnitude  square 
in  step  (iv)  is  detailed  in  section  4.2.  Determination  of  the  system  impulse  functions  and  phases 
required  in  step  (v)  are  given  respectively  in  sections  4.3  and  4.4.  Computation  of  the  emp 


response  of  the  same  system  is  derived  in  section  4.5.  Related  numerical  results  are  presented  in 
figures  4  through  52. 

The  final  fi^equency  and  time  responses  obtained  from  this  method  may  then  be  analyzed  and 
compared  with  the  measured  data,  if  available. 


2.  Approximation  of  Resonances  by  Second-Order  Transfer  Functions 

The  second-order  transfer  function  used  to  approximate  each  cw  resonance  displayed  on 
measured  magnitude  data  takes  either  of  the  following  two  forms: 

^ci'">  -  ^^ '  (la) 

s    ^  as  ^  b 

or 

jj  ..  A{s  ^  c) 

s    ^  as  *  b 

where  s  is  the  complex  frequency,  a,  b  and  A  are  all  positive  and  real  numbers  with  b  >  a^ll  [1], 
and  c  is  a  real  number  (which  may  be  positive,  0,  or  negative). 

For  the  transfer  function  in  eq  (la),  the  parameters  a,  b  and  A  are  determined  by  [1] 

b'-   <^l-        \     ''    .  (2a) 
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a'  -  2{b   -  G)o)  ,  (2b) 

and 

A  =  {(D^  -  (Oj) ,  (2c) 

where  Wq  is  the  radian  resonant  frequency,  Wj  and  co,  are  the  half-power  frequencies  with 
Wi  <  Wo  <  0)2  and  \H^(J^)\  =  WJJ^^ )\  =  \HJJ^^) ^\l ^  ,  and  |//^(/Wo  )|  is  the  maximum 
magnitude  at  the  resonant  frequency. 

All  of  the  four  required  numbers,  Wg ,  Oi  ,  0)2 ,  and  \H^(Jwq  )|,  can  be  read  directly  from  the 
measured  cw  magnitude  data.  In  this  case,  (Jj  and  Wj  are  not  independent,  but  are  related  by 

Wj     +    (x>2   =     2(x)q  (3) 

Thus,  it  is  unimportant  if  either  coi  or  0)2  is  not  available  in  the  measured  curve  since  it  can  be 
determined  from  eq  (3).  The  dimensions  of  a  and  b  in  this  case  are  respectively  s'^  and  s'^.  The 
dimension  of  ^  depends  on  |//^(/a))|. 


For  the  transfer  function  eq  (lb),  the  parameters  a,  b,  c,  and  A  are  determined  by  [1] 

b^  =  2g)J  -   cOjWj  ,  (4a) 


2  ,2  _2  .2 


a^  =  Qj  +   G)2  -  4a)o  +  26  ,  (4b) 


2       2s 


(Wq    -     "l<»^2) 

^     =  —, -, '  (4c) 

(g)i  .  (O2  -  2a)o) 

and 

A'  =(2co'o  -2Z)  .a^)|//,0«J|^  .  (4d) 

Again,  the  required  input  information,  cOq  ,  coi  ,  Wj,  and  \HJJix>q  )|,  can  be  read  directly  from  the 
given  curve.  The  dimensions  of  a,  b  and  c  now  are  respectively  s'\  s'^  and  s"\  The  dimension  of 
A  depends  on  \Hi,(joiQ  )|. 

In  eq  (4c),  we  require  that  [1] 

2  2  2  2 

0)^  +  0)  2  >  2  G)(,      and       o)  q  >  co  ^  0)2   ■  (5) 

If  either  of  the  conditions  in  eq  (5)  is  not  satisfied,  we  use  the  transfer  function  in  eq  (la). 

Once  the  parameters  a,  b,  c  [under  condition  eq  (5)],  and  A  are  determined,  the  magnitude  square 
of  the  component  transfer  function  that  approximately  represents  the  resonant  frequency  under 
consideration  becomes 

\H^(J^)\'  - 4— ,  (6a) 

or 

I   IT-/    •      M2  A^(^x>^     +    C^) 

CO    -  {lb  -  a   )(x)    +  b 

The  approximation  as  presented  is  good  only  for  the  important  resonant  region  being  studied  and 
is  not  meant  to  match  the  given  measured  curve  everywhere. 

If  there  are  A/^  distinct  resonant  frequencies  in  the  original  cw  data,  the  final  approximate 
magnitude  square  for  the  overall  system  is  just  the  sum  of  all  the  component  magnitude  squares, 
from  which  the  transfer  function  in  complex  frequency  H(s)  is  then  extracted  [2].  It  will  be 
shown  later  that  there  may  be  multiple  solutions  for  H(s).  The  corresponding  solutions  for 
impulse  transfer  fijnctions  h(t)  can  then  be  determined  in  terms  of  damped  sine  and  cosine 
functions. 


3.  Frequency  Transformation 

To  avoid  manipulations  with  large  numbers  such  as  the  frequency  in  megahertz  (MHz),  we  use 
the  normalized  radian  frequency  by  omitting  the  factor  10^.  We  herein  designate  o)  as  the 
normalized  radian  frequency  and  to'  as  the  actual  radian  frequency  such  that 

G)'  =  5(0         with         B  =  \0^  .  (7) 

Substituting  eq  (7)  in  eqs  (6a)  and  (6b),  we  obtain  respectively  the  component  magnitude  square 
in  actual  radian  frequency. 


2  d4 


\G^U^') 


i\\2 


A^B 


and 


(0''  -  B\2b  -  a')(o''  .  b^B* 

£^ 

CO''  -  (lb'  -  a'')(o''  .  b'^ 


0)'*  -B^{2b  -  a')(o''  .  b^B' 
G)''  -  (lb'  -  a''')oi'''  .  *'' 


(8a) 


(8b) 


?2     AI>2 


From  eq  (8a),  we  see  that  the  parameters  a\  b\  and  A'  will  become  respectively  aB,  bB  ,  AB 
when  we  convert  from  normalized  frequency  to  actual  frequency.  For  the  transfer  function  in  eq 
(lb),  the  parameters  a\  b\  c\  and  A'  in  eq  (8b)  will  become  aB,  bB^,  cB,  and  AB  for  the  same 
conversion  in  frequency. 


4.  Assessment  of  Torpedo  Measurement  Data 

The  data  analyzed  in  this  report  came  from  a  test  made  on  a  special  torpedo  by  Naval  Surface 
Warfare  Center  (NSWC).  The  test  was  instrumented  with  current  probes  in  two  locations, 
arbitrarily  called  C26  and  C29,  where  EMI  problems  are  suspected.  The  external  cw  source  used 
in  the  test  for  C26  is  a  magnetic  field  in  amperes/meter  (input)  as  shown  in  figure  1 .  The 
measured  probe  magnitude  response  in  amperes  (output)  is  given  in  figure  2.  Thus,  this  input- 
output  configuration  constitutes  "the  system"  as  far  as  our  particular  application  is  concerned. 
Both  figures  1  and  2,  were  originally  expressed  in  megahertz.  For  our  analysis,  we  use  the 
normalized  frequency  by  omitting  the  factor  of  10^,  as  discussed  in  section  3.  The  transfer 
function  magnitude  in  meters  is  obtained  by  taking  the  ratio  of  the  response  data  (fig.  2)  to  the 
source  data  (fig.  1),  as  shown  in  figure  3  (solid  curve),  where  the  normalized  frequency  is  used. 
Also,  we  concentrate  only  on  the  portion  of  data  up  to  100  Hz  for  the  normalized  frequency  or 
100  MHz  in  real  frequency.  Similar  analysis  for  C29  is  not  included  in  this  report. 


4.1.  Analytical  Expression  for  Magnitude  Squares 

From  the  measured  data  derived  in  figure  3  we  see  many  resonant  frequencies  where  the 
measured  magnitudes  reach  their  maxima.  For  demonstration  and  considering  the  storage 
limitation  of  our  software,  as  a  first  step,  we  identify  six  resonant  frequencies  properly  numbered 
in  the  figure  3.  They  are  :  fjo  =  5.593,  fjo  =  19.609,  fjo  =  53.262,  f4o  =  56.848,  fjo  =  59.373,  and 
fgo  =  72.192  with  their  respective  maximum  values  of  0.1 151  m,  0.0862  m,  0.3773  m,  0.2709  m, 
0.4671m,  and  0.1277  m. 

For  the  first  resonant  fi^equency  at  fjo  (cojo  =  35. 143 1),  we  observe  the  half-power  frequencies 
fji  =5.077  (g)ii  =  31.9016)  and  fi2  =5.935  (<0i2  =37.2876).  Since  one  of  the  conditions  in  eq 
(5)  is  not  satisfied  (wn^  +  Wij^  <  2a)io^ ),  we  use  the  transfer  function  in  eq  (la)  and  obtain  a^ 
=  5.2870,  Z>i  =  1.2490(10^),  mdA^  =  21.4460.  Thus,  the  first  component  magnitude  square 
becomes 

WSJ^^)?  -  '—  ,  (9a) 

where 

Z)j(o)2)  =  w'  -  2.4701(10^)0)2  .  1.5600(10^)  .  (9b) 

The  magnitude  square  in  eq  (9a)  is  used  to  approximate  the  first  resonant  region  in  figure  3. 

For  the  second  resonant  fi-equency  at  fjo  =  19.609  (wjo  =  123.2070),  we  read  fji  =  17.784  (oji 
=  1 1 1.7402),  and  fjj  =  20.480  (0)22  =  128.6796).  Again,  because  of  0)21^  +  022^  <  20)20^ ,  we  use 
the  transfer  function  in  eq  (la)  to  obtain  aj  =  16.4905,  b^  =  1.5316(10^),  A^  =  1.7545(10^),  and 
the  second  component  magnitude  square 

\H0(^)\^  =  —  ,  (10a) 

i),(a>2) 

where 

Dl^i^)  =   (^'  -  3.0360(10^)0)2  .  2.3458(10')  .  (10b) 

For  the  third  resonant  frequency  at  fjo  =  53.262  (0)30  =  334.6550),  we  have  fji  =  52.604  (0)31 
=  330.5207),  and  f32  =  55.124  (0)32  =  346.3543).  Since  the  second  condition  in  eq  (5)  is  not 
satisfied,  that  is  0)30^  <  o)3j  0)32 ,  we  also  use  the  transfer  function  in  eq  (la)  to  obtain  a^ 
=  16.0080,  Z>3  =  1.1212(10'),  ^3  =  2.0218(10^),  and  the  third  component  magnitude  square 

I^3(/'")r=   ^V'  Ola) 

Z)3(0)2) 

where 

Z)3(o)2)  =   0)'  -  22399(lCP)o)2  .  1.2571(10'°)  .  (lib) 


For  the  fourth  resonant  frequency  at  f4o  =  56.848  (040  =  357. 1865),  where  f4i  does  not  exist  and 
f42  =  58.324  (042  =  366.4605),  we  also  use  eq  (la)  to  obtain  a^  =  18.7740,  b^  =  1.2776(10^), 
^4  =  1.8172(10^),  and  the  fourth  component  magnitude  square 

m<^)\'-  -^,  (12a) 

where 

i)/a)2)  =  co"  -  2.5516(10^)0)'  .  1.6322(10'°)  .  (12b) 

Similarly,  for  the  fifth  and  sixth  resonance  at  fjo  =  59.373  and  f^  =  72.192,  the  transfer  function 
in  eq  (la)  also  applies.  We  obtain  a^  =  9.6223,  b^  =  1.3921(10^),  A^  =  1.6769(10^);  a^  =  19.2431, 
b^  =  2.0593(10^),  and  ^6  =  1.1149(10^);  which  yield  respectively 

|//3(/G))|'=  ^i— ,  (13a) 


with 


and 


Z)/a)')  =  (0^  -  27833(lCP)a)'  .  1.9380(10'°),  (I3b) 


Al 


with 


\H,(j^W-   -^-,  (14a) 


D/o)')  =  G)'  -  4.1150(10?)g)2  .  4.2409(10'°)  .  (Mb) 


Adding  all  the  component  magnitude  squares  in  eqs  (9a)  through  (14a),  we  have  the  overall 
system  magnitude  square, 

6 

\H(jcor-  y\KU^)\'  -   1.1473(1C7)^^^  ,  (15) 


where 


with  Z),  (g)'),  i  =  1,  2,  ...,  6  given  previously  in  eqs  (9b)  to  (14b),  and 


iV(a)2)  =  0)'°  -  9.3575(10^)0)''  .  3.6725(10")o)''  -  7.8083  (10")o)" 

.  9.6766(10"  )o)''  -  6.9726(10'')o)'°  .  2.7329(10'')o)'  -  5.0434(lCP')o)' 
.  4.0175(10'')o)'  -  8.0753(10*^)0)'  .  6.0670(10*') 


=  [o'  -  2.4904(10')(o'  +   2.1141(10')]   [w'  -  3.0973(10^"  -  2.7749(10')] 

♦  [o)'  -  23947(10')w'  .  1.4446(1^°)]  [ui' -  2.6894(10')co' .    1.8147(lC/° )]      (17) 

*  [o*  -  3.9387(10^)0)'  .  3.9452(l(y°)]  . 

The  final  maximum  value  of\H(joi)\  in  eq  (15)  at  Wjo  is  0.4921,  somewhat  larger  than  the  actual 
maximum  value  of  0.4671,  resulting  in  an  amplification  factor  of  1 .0535.  This  increase  is  due  to 
the  contributions,  although  small,  from  the  other  component  magnitude  squares.  To  enforce 
|/f(/a)5o)|  =  0.4671,  we  simply  divide  the  final  magnitude  square  in  eq  (15)  by  (1.0535)'  to  yield  a 
corrected  system  magnitude  square 

|///»|'  =  1.0336(10^)^^,  (18) 

with  D(iii^)  and  iV(a)')  given  in  eqs  (16)  and  (17). 

The  square  root  of  eq  (18)  is  plotted  as  the  broken  curve  in  figure  3  for  comparison.  As  far  as 
the  major  regions  (resonant  frequencies)  in  figure  3  are  concerned,  the  approximation  is  very 
good  except  that  the  approximated  value  at  f4o  is  much  larger  than  the  actual  value. 

According  to  the  method  known  in  classical  network  theory  [2],  we  can  extract  the  transfer 
function  in  complex  frequency  s  from  the  magnitude  square  in  eq  (18)  by  substituting  s  =ji^. 
Thus, 

H{s)H{-s)  -  \H(ji^)\^        {withs-ji^) 


1.0336(10')^t^ 


where 


and 


with 


N{-s^)  -  N^i^)Nl-)N^(^)Nf)N^{.)Nl-)N^{^)Nl-)N^(.)Nl-), 
D{-s')  -  D^{.)D^i-)D2(^)Df)Dl^)D,i-)D^(.)D^i-)Dl.)Dl-)Dl.)D^i-) 


N^(±)  =  s^  ±  20.43435  .  1.4540(10^) 
Nl±)  =  s^  ±  48.40785  .  1.6658(10^ 
NJi±)  =  5'  ±  30.13245  .  1.2019(10') 
A^/±)  =  5'  ±  21.92905  .  1.3471(10') 
Nji±)  =  5'  ±  58.15255  +  1.9862(10') 


(19) 

(20) 
(21) 


(22) 


and 


Z).(±)  =  5'  ±  as  +  b. 


7=   1,2,. ..,6  . 


(23) 


4. 2.  Extraction  of  Transfer  Functions 

To  determine  a  realizable  transfer  function  H{s)  from  eq  (19),  we  have  to  assign  all  D,{+)  to  the 
denominator  ofH(s)  in  order  to  yield  a  stable  system.  Then  all  Z),(-)  automatically  belong  to  the 
denominator  ofH{-s).  As  far  as  the  numerator  of  H(s)  is  concerned,  we  have  multiple  choices. 
If  all  of  the  Ni(+)  are  chosen  as  the  numerator  of  H(s),  Nj(-)  will  be  in  the  numerator  of  H{-s). 
In  this  case,  H(s)  has  no  zeros  in  the  right-half  s-plane.  The  system  transfer  function  is  said  to  be 
at  minimum  phase.  Other  choices  will  produce  one  or  more  zeros  in  the  right-half  s-plane, 
resulting  in  nonminimum-phase  transfer  functions.  There  are  a  total  of  3 1  (from  2^  - 1) 
nonminimum-phase  transfer  functions  for  the  case  under  study.  The  coeflQcient  in  eq  (19), 
1.0334(10')  gives  a  factor  of  3.2146(10^)  each  to  His)  and  H{-s). 


For  the  minimum-phase  case,  we  obtain 

//i5)  =  3.2146(10' )-ii 


n^(-) 


n^,o 


i-l 


where 


yKs  .  G 

L    Di.)      ' 

Fr  = 

0.4305, 

Gi  =  4.5011, 

F2  = 

0.8814, 

G2  =  46.740, 

Fs  = 

2.5998, 

G3=  1196.6, 

F,= 

-0.1656, 

G4  =  950.33, 

Fs  = 

-2.0624, 

G5  =  926.21, 

Fe- 

-1.6837, 

G6  =  93.272. 

(24) 


(25) 
and 


If  we  replace  one  or  more  Ni(+)  in  eq  (24)  by  its  counterpart  N,(-),  the  resuhing  transfer  function 
wiU  be  at  nonminimum  phase,  and  the  corresponding  expansion  coefiQcients  F^  and  G,  will  be 
different  from  those  in  eq  (25).  The  values  of  F,  and  G,  for  the  3 1  nonminimum-phase  cases  can 
be  computed  and  printed  out  by  the  software,  when  desired. 


4.3.  Impulse  Response 

The  impulse  response  for  the  minimum-phase  case  in  eq  (24)  is  then 

K(^)  -  )\,cosi^,t)  .  ^,sin(p,.0]^""''  ,  (26) 


with 


a,  = 

2.6435, 

P,  =  35.243, 

c,  =  0.4305, 

d,= 

0.0954, 

a?  = 

8.2443, 

P2=  123.48, 

C2=  0.8814, 

d,- 

0.3197, 

"3  = 

8.0012, 

p3  =  334.75, 

C3=  2.5998, 

d,= 

3.5126, 

"4  = 

9.3860, 

P4  =  357.31, 

C4  = -0.1656, 

d,= 

2.6640, 

a,= 

4.8106, 

P5  =  373.08, 

C5  =  -2.0624, 

ds  = 

2.5092, 

and 

"6  = 

9.6219, 

P6  =  453.70, 

C6  = -1.6837, 

de  = 

0.2129. 

(27) 


Equation  (26)  is  plotted  in  figure  4  marked  with  '11111',  because  five  N^  (+)  are  included  in 
the  numerator  of  eq  (24)  to  represent  the  minimum-phase  case.    The  time  scale,  in  seconds,  used 
here  is  based  on  the  normalized  frequency.  The  corresponding  time  scale  for  the  actual 
frequency  will  be  microseconds.  Nonminimum  cases  are  shown  in  figures  4  to  1 1  for 
comparison.  The  coefficients  c,  and  d^  for  the  nonminimum-phase  cases  will  all  be  different  from 
those  in  eq  (27)  although  a,  and  P,  will  remain  unchanged.  The  curve  marked  with  '1  1  - 1  - 1  1' 
is  interpreted  to  indicate  that  A^i(+)  A'2(+)  A3(-)  A^4(-)  ^si"^)  is  chosen  for  the  numerator  of  the 
transfer  function. 


4.4.  Phase  Information 

Since  we  define,  in  this  report,  the  phase  of //(/to)  as 

Hgco)  -  |//Oa>)k  ^«(">  , 
we  have  for  the  minimum-phase  case  in  eq  (24) 

0^(0))    =    0/0))    -    0„(G))    , 

where  0/w),  due  to  the  denominator  factors  in  H^  (/co),  is  given  by 


(28) 


(29) 


0/0))  =    >  tan-^ 


6 


»i 


a  0) 


0) 


(30) 


and  0„(o)),  due  to  the  numerator  factors  in  //,„(/o)),  from  eq  (22),  is  given  by 


e„(a))  =  tan"' 


20.4343  0) 


1.4540(10^)   -  0)2 


tan' 


48.4078  0) 


1.6658(10'*)  -  0)2 


tan 


30. 1324  0) 


1.2019(10^)   -  0)2 


+  tan" 


21.9290  0) 


1.3471  (10^)   -  0)2 


tan" 


58.1525  0) 


1.9862(10^)   -  0)2 


(31) 


For  the  nonminimum-phase  cases,  one  or  all  of  the  terms  in  eq  (3 1)  will  change  sign  to  make  the 
total  phase  larger  than  that  for  the  minimum-phase  case.  Phases  for  the  32  cases  (1  minimum 
phase  and  31  nonminimum  phases)  are  presented  in  figures  12  through  19. 


4.5.  EMP  Response 

Once  the  impulse  response  h{t)  is  obtained  in  subsection  4.3,  the  response  r^^Jif)  at  C26  due  to 
an  emp  source  function  emp{i)  can  be  derived  analytically  by  performing  convolution  integrals. 
That  is, 

'•.JO   =  ^^P^^)  *  ^0  ,  (32) 

where  h{t)  represents  either  minimum-phase  or  nonminimum-phase  case. 

The  actual  emp  source  function  used  in  the  test  is  given  by 


emp{t')  =  K 


g-4(10*)r'   _    g  475(10«)r' 


where  t  is  the  real  time  before  frequency  normalization,  and  the  constant  factor  K  is 
J,=    52500^   139.2573. 

377 
The  Laplace  transform  of  eq  (33)  is 


(33) 


EMP{s')  =  K 


1 


1 


s'  .  4(10^)       s'  .  475(10') 


(34) 


and  its  spectrum  is  given  by 


EMPiJoi')  =  K 


1 


1 


4(10')  .70'        475(10')  .  >' 


(35) 


We  now  apply  frequency  normalization  to  convert  eq  (35)  into 


EMPHio)  -K(m') 


1 


1 


4  +  j(x>        475  +  y'o) 


which  yields 


emp^  (0  =  K{\Cr^)[e"  -  e''"] . 


(36) 


(37) 


Comparing  eqs  (33)  and  (37),  we  see  that  the  net  effects  on  the  time  excitation  function  from  the 
frequency  normalization  are  a  dropping  of  the  (10')  factor  from  the  exponent  and  a  multiplying 
factor  of  (10"')  to  the  coefficient. 

Now,  the  frequency  response  to  the  normalized  source  becomes,  from  eq  (36), 


10 


l^(/")i  =  \Hin^)\- ,  (38) 

v/(a)'  .  4')(a)'  .  475' ) 

which  is  plotted  in  figure  20.  Since  there  is  only  one  |//c(/w)|  for  the  minimum-phase  and  non- 
minimum-phase  cases,  there  is  only  one  fi-equency  response  as  shown  in  figure  20. 

The  time  response  to  the  normalized  source  becomes 

r^J)  -  Kim')^''  -  e-'''^].  Kt)  (39) 

Since  the  source  fianction  is  in  the  form  of  e'^'  and  the  typical  term  in  h{f)  is 
]g.  cos(  P  .  0   +  d.  sin(  p.  0  ]  ^  "'  ,  we  only  need  to  deal  with  a  general  formula 

e*'  *  [c.cos(p.O  +  J.sin(P7)  V"''.  One  term  is 


t 


e-''  *  cos(p.O  e    '    =      e-"^'-'^  cos(p.x)  e  '    dx 


[{p  -   a^)cos(p,0   .  p,sm(p7)]e-"''  -  {p  -  a,>-'' 


{p  -  cc,T  ^  P? 
and  the  other  term  is 


e-''  *   sin(p.O  e    '    =      e^^'-^^  sin(p.T)  e    '    dx 


[-P,cos(P,0  .  {p   -  a^)sin(P,0]^""''  ^  P,^"' 


ip  -  ay  .  p,.^ 

Substituting  eqs  (40a)  and  (40b)  into  eq  (39),  we  obtain 


i  )  l/?2.  'e ""''  +  «.' 


where 

,       jc^p^  -  g,)  -  ^,pjcos(p,0  .   [c  p,  .  ^,(/>,   g  )]sm(p,0 

(/'i-   a)'  .  P^ 
and 

11 


(40a) 


(40b) 


(41) 


(42a) 


-c(p    -   a)  +  d^ 
«/ =   — '■ '■ '—!-  {independent of  t)  ,  /'42k\ 

and  m"  and  /?,"  have  the  same  expressions  as  w/  and  n-  but  withpi  replaced  by/?2 ,  where/?,  =  4 
andpj'^^TS. 

The  time  responses  for  all  cases  are  presented  in  figures  21  to  52,  which  can  be  compared  with 
the  measured  emp  response  as  shown  in  figure  53.  The  emp  response  r^^Jif)  for  the  minimum- 
phase  case  as  shown  in  figure  2 1  is  almost  always  nonnegative  because  the  term  in  eq  (42b)  is 
very  dominant  in  this  case.  Similar  solutions  for  the  nonminimum-phase  cases  can  be  both 
positive  and  negative.  Some  of  them  have  significant  negative  parts  such  as  that  shown  in  figure 
38.  Therefore,  it  appears  that  the  specific  input-output  configuration  of  the  torpedo  under  study 
constitutes  more  likely  a  nonminimum-phase  system. 


5.  Conclusions 

We  have  applied  the  newly  developed  theory  for  predicting  the  impulse  characteristics  of  an 
unknown  linear  system  based  only  on  available  cw  magnitude  data  to  a  particular  system 
configuration.  By  comparing  the  predicted  impulse  responses  to  the  measured  emp  data,  we 
conclude  that  the  system  under  study  is  more  likely  at  nonminimum  phase.  This  same  procedure 
can  be  applied  to  the  analysis  of  other  linear  systems  and  to  assess  the  meaning  of  the  measured 
data. 


We  thank  the  Naval  Surface  Warfare  Center  (NSWC)  and  Naval  Sea  Systems  Command 
(NAVSEA)  for  their  partial  support  of  this  study.  We  also  appreciate  the  measured  data  supplied 
by  John  Bean  and  Rich  Porter  of  NSWC. 


6.  References 

[1]        Ma,  M.T.;  Adams,  J.W.  Characteristics  of  unknown  linear  systems  deduced  from 

measured  cw  magnitude.  J.  Res.  Natl.  Inst.  Stand.  Technol.,  98:297-319;  May- June 
1993. 

[2]        Balabanian,  N.  Network  Synthesis,  New  Jersey:  Prentice-Hall;  1958. 


12 


20 


o 

(D 
C 
O) 
CO 


15 


E 
E^ 

■^       10 


5  — 


20  40  60 

Frequency  (MHz) 


80 


100 


Figure  1. 


The  cw  external  excitation  (magnetic  field)  used  in  the  experiment. 


13 


1.6 


1.2  — 


< 
E 

-£       0.8 
(D 

o 


0.4 


0.0 


20  40  60 

Frequency  (MHz) 


80 


100 


Figure  2.  Measured  C26  magnitude  response  to  the  cw  excitation  shown  in  Figure  1 . 
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Figure  3.  Magnitude  of  the  system  transfer  function,  both  measured  and 

approximated  based  on  eq  (18). 
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Figure  4.  Possible  impulse  responses  of  the  system  under  study;  minimum  phase  case 

and  three  nonminimum  phase  cases. 


16 


(U 

■o 

•^ 
Q. 

E 
< 


Curve  5  = 

1 

1 

-1 

1 

1 

Curve  6  = 

1 

1 

-1 

1 

-1 

Curve  7  = 

1 

1 

-1 

-1 

1 

Curve  8  = 

1 

1 

-1 

-1 

-1 

Curve  7 


0) 

TJ 

3 
+^ 

Q. 

E 
< 


Curves 


Figure  5.  Possible  impulse  responses  of  the  system  under  study;  four  nonminimum 
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Figure  7.  Possible  impulse  responses  of  the  system  under  study,  four  nonminimum 

phase  cases. 
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Figure  8.  Possible  impulse  responses  of  the  system  under  study,  four  nonminimum 

phase  cases. 
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Figure  9.  Possible  impulse  responses  of  the  system  under  study,  four  nonminimum 

phase  cases. 
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Figure  10.        Possible  impulse  responses  of  the  system  under  study,  four  nonminimum 
phase  cases. 
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Figure  1 1 .        Possible  impulse  responses  of  the  system  under  study,  four  nonminimum 
phase  cases. 
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Figure  12.        Possible  phases  of  the  system  under  study  (C26)  for  one  minimum  phase 
case  and  three  nonminimum  phase  cases. 
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Figure  13.        Possible  phases  of  the  system  under  study  (C26)  of  four  nonminimum 
phase  cases. 
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Figure  14.        Possible  phases  of  the  system  under  study  (C26)  of  four  nonminimum 
phase  cases. 
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Figure  15.        Possible  phases  of  the  system  under  study  (C26)  of  four  nonminimum 
phase  cases. 
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Figure  16.        Possible  phases  of  the  system  under  study  (C26)  of  four  nonminimum 
phase  cases. 
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Figure  17.        Possible  phases  of  the  system  under  study  (C26)  of  four  nonminimum 
phase  cases. 
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Figure  18.        Possible  phases  of  the  system  under  study  (C26)  of  four  nonminimum 
phase  cases. 
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Figure  19.        Possible  phases  of  the  system  under  study  (C26)  of  four  nonminimum 
phase  cases. 
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Figure  20.        Magnitude  response,  as  a  function  of  normalized  frequency,  of  the  system 
to  the  impulse  excitation  given  in  eq  (36). 
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Figure  21.        Time  response  of  the  system  to  the  impulse  excitation  given  in  eq  (37), 
for  the  minimum  phase  case. 
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Figure  22.        Time  response  of  the  system  to  the  impulse  excitation  given  in  eq  (37), 
for  a  nonminimum  phase  case. 
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Figure  23.        Time  response  of  the  system  to  the  impulse  excitation  given  in  eq  (37), 
for  a  nonminimum  phase  case. 
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Figure  24.        Time  response  of  the  system  to  the  impulse  excitation  given  in  eq  (37), 
for  a  nonminimum  phase  case. 
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Figure  25.        Time  response  of  the  system  to  the  impulse  excitation  given  in  eq  (37), 
for  a  nonminimum  phase  case. 
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Figure  26.        Time  response  of  the  system  to  the  impulse  excitation  given  in  eq  (37), 
for  a  nonminimum  phase  case. 
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Figure  27.        Time  response  of  the  system  to  the  impulse  excitation  given  in  eq  (37), 
for  a  nonminimum  phase  case. 
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Figure  28.        Time  response  of  the  system  to  the  impulse  excitation  given  in  eq  (37), 
for  a  nonminimum  phase  case. 
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Figure  29.        Time  response  of  the  system  to  the  impulse  excitation  given  in  eq  (37), 
for  a  nonminimum  phase  case. 
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Figure  30.        Time  response  of  the  system  to  the  impulse  excitation  given  in  eq  (37), 
for  a  nonminimum  phase  case. 
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Figure  3 1 .        Time  response  of  the  system  to  the  impulse  excitation  given  in  eq  (37), 
for  a  nonminimum  phase  case. 
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Figure  32.        Time  response  of  the  system  to  the  impulse  excitation  given  in  eq  (37), 
for  a  nonminimum  phase  case. 
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Figure  33.        Time  response  of  the  system  to  the  impulse  excitation  given  in  eq  (37), 
for  a  nonminimum  phase  case. 
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Figure  34.        Time  response  of  the  system  to  the  impulse  excitation  given  in  eq  (37), 
for  a  nonminimum  phase  case. 
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Figure  35.        Time  response  of  the  system  to  the  impulse  excitation  given  in  eq  (37), 
for  a  nonminimum  phase  case. 
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Figure  36.        Time  response  of  the  system  to  the  impulse  excitation  given  in  eq  (37), 
for  a  nonminimum  phase  case. 
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Figure  37.        Time  response  of  the  system  to  the  impulse  excitation  given  in  eq  (37), 
for  a  nonminimum  phase  case. 
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Figure  38.        Time  response  of  the  system  to  the  impulse  excitation  given  in  eq  (37), 
for  a  nonminimum  phase  case. 
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Figure  39.        Time  response  of  the  system  to  the  impulse  excitation  given  in  eq  (37), 
for  a  nonminimum  phase  case. 
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Figure  40.        Time  response  of  the  system  to  the  impulse  excitation  given  in  eq  (37), 
for  a  nonminimum  phase  case. 
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Figure  41 .        Time  response  of  the  system  to  the  impulse  excitation  given  in  eq  (37), 
for  a  nonminimum  phase  case. 
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Figure  42.        Time  response  of  the  system  to  the  impulse  excitation  given  in  eq  (37), 
for  a  nonminimum  phase  case. 
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Figure  43.        Time  response  of  the  system  to  the  impulse  excitation  given  in  eq  (37), 
for  a  nonminimum  phase  case. 


55 


curve  =  24 


(-1  1  -1  -1  -1  ) 


a  = 


2.64345 
8.24428 
8.00117 
9.38598 
4.81063 
9.62187 


P  = 


35.2426 
123.482 
334.749 
357.312 
373.079 
453.699 


c  = 


-0.711 

0.96383 

-9.52379 

-1.50512 

7.52899 

3.24709 


d  = 


0.23855 

0.17798 

-3.04781 

14.3331 

-1.59108 

-0.74954 


3 


Time  (s) 


Figure  44.        Time  response  of  the  system  to  the  impulse  excitation  given  in  eq  (37), 
for  a  nonminimum  phase  case. 
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Figure  45.        Time  response  of  the  system  to  the  impulse  excitation  given  in  eq  (37), 
for  a  nonminimum  phase  case. 
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Figure  46.        Time  response  of  the  system  to  the  impulse  excitation  given  in  eq  (37), 
for  a  nonminimum  phase  case. 
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Figure  47.        Time  response  of  the  system  to  the  impulse  excitation  given  in  eq  (37), 
for  a  nonminimum  phase  case. 
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Figure  48.        Time  response  of  the  system  to  the  impulse  excitation  given  in  eq  (37), 
for  a  nonminimum  phase  case. 
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Figure  49.        Time  response  of  the  system  to  the  impulse  excitation  given  in  eq  (37), 
for  a  nonminimum  phase  case. 
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Figure  50.        Time  response  of  the  system  to  the  impulse  excitation  given  in  eq  (37), 
for  a  nonminimum  phase  case. 
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Figure  5 1 .        Time  response  of  the  system  to  the  impulse  excitation  given  in  eq  (37), 
for  a  nonminimum  phase  case. 
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Figure  52.        Time  response  of  the  system  to  the  impulse  excitation  given  in  eq  (37), 
for  a  nonminimum  phase  case. 
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Figure  53.        Measured  time  response  of  C26,  after  normalization,  to  the  same  impulse 
excitation. 
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Technical  Publications 

Periodical 

Journal  of  Research  of  the  National  Institute  of  Standards  and  Technology— Reports  NIST  research 
and  development  in  those  disciplines  of  the  physical  and  engineering  sciences  in  which  the  Institute  is 
active.   These  include  physics,  chemistry,  engineering,  mathematics,  and  computer  sciences.    Papers 
cover  a  broad  range  of  subjects,  with  major  emphasis  on  measurement  methodology  and  the  basic 
technology  underlying  standardization.    Also  included  from  time  to  time  are  survey  articles  on  topics 
closely  related  to  the  Institute's  technical  and  scientific  programs.    Issued  six  times  a  year. 

Nonperiodicals 

Monographs— Major  contributions  to  the  technical  literature  on  various  subjects  related  to  the  Institute's 
scientific  and  technical  activities. 

Handbooks— Recommended  codes  of  engineering  and  industrial  practice  (including  safety  codes) 
developed  in  cooperation  with  interested  industries,  professional  organizations,  and  regulatory  bodies. 

Special  Publications— Include  proceedings  of  conferences  sponsored  by  NIST,  NIST  annual  reports,  and 
other  special  publications  appropriate  to  this  grouping  such  as  wall  charts,  pocket  cards,  and 
bibliographies. 

Applied  Mathematics  Series— Mathematical  tables,  manuals,  and  studies  of  special  interest  to  physicists, 
engineers,  chemists,  biologists,  mathematicians,  computer  programmers,  and  others  engaged  in  scientific 
and  technical  work. 

National  Standard  Reference  Data  Series— Provides  quantitative  data  on  the  physical  and  chemical 
properties  of  materials,  compiled  from  the  world's  literature  and  critically  evaluated.    Developed  under  a 
worldwide  program  coordinated  by  NIST  under  the  authority  of  the  National  Standard  Data  Act  (Public 
Law  90-396).    NOTE:   The  Journal  of  Physical  and  Chemical  Reference  Data  (JPCRD)  is  published  bi- 
monthly for  NIST  by  the  American  Chemical  Society  (ACS)  and  the  American  Institute  of  Physics  (AIP). 
Subscriptions,  reprints,  and  supplements  are  available  from  ACS,  1155  Sixteenth  St.,  NW,  Washington, 
DC   20056. 

Building  Science  Series— Disseminates  technical  information  developed  at  the  Institute  on  building 
materials,  components,  systems,  and  whole  structures.   The  series  presents  research  results,  test 
methods,  and  performance  criteria  related  to  the  structural  and  environmental  functions  and  the  durability 
and  safety  characteristics  of  building  elements  and  systems. 

Technical  Notes— Studies  or  reports  which  are  complete  in  themselves  but  restrictive  in  their  treatment 
of  a  subject.    Analogous  to  monographs  but  not  so  comprehensive  in  scope  or  definitive  in  treatment  of 
the  subject  area.    Often  serve  as  a  vehicle  for  final  reports  of  work  performed  at  NIST  under  the 
sponsorship  of  other  government  agencies. 

Voluntary  Product  Standards— Developed  under  procedures  published  by  the  Department  of  Commerce 
in  Part  10,  Title  15,  of  the  Code  of  Federal  Regulations.   The  standards  establish  nationally  recognized 
requirements  for  products,  and  provide  all  concerned  interests  with  a  basis  for  common  understanding  of 
the  characteristics  of  the  products.    NIST  administers  this  program  in  support  of  the  efforts  of  private- 
sector  standardizing  organizations. 

Consumer  Information  Series— Practical  information,  based  on  NIST  research  and  experience,  covering 
areas  of  interest  to  the  consumer.    Easily  understandable  language  and  illustrations  provide  useful 
background  knowledge  for  shopping  in  today's  technological  marketplace. 

Order  the  above  NIST  publications  from:  Superintendent  of  Documents,  Government  Printing  Office, 
Washington.  DC  20402. 

Order  the  following  NIST  publications— FIPS  and  NISTIRs—from  the  National  Technical  Information 
Service,  Springfield,  VA  22161. 

Federal  Information  Processing  Standards  Publications  (FIPS  PUB)— Publications  in  this  series 
collectively  constitute  the  Federal  Information  Processing  Standards  Register.   The  Register  serves  as  the 
official  source  of  information  in  the  Federal  Government  regarding  standards  issued  by  NIST  pursuant  to 
the  Federal  Property  and  Administrative  Services  Act  of  1949  as  amended.  Public  Law  89-306  (79  Stat. 
1127),  and  as  implemented  by  Executive  Order  11717  (38  FR  12315,  dated  May  11,  1973)  and  Part  6  of 
Title  15  CFR  (Code  of  Federal  Regulations). 

NIST  Interagency  Reports  (NISTIR)— A  special  series  of  interim  or  final  reports  on  work  performed  by 
NIST  for  outside  sponsors  (both  government  and  non-government).    In  general,  initial  distribution  is 
handled  by  the  sponsor;  public  distribution  is  by  the  National  Technical  Information  Service,  Springfield, 
VA  22161,  in  paper  copy  or  microfiche  form. 


